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Abstract. The Poincare gauge theory of gravity has a metric compatible connection with independent 
dynamics that is reflected in the torsion and curvature. The theory allows two good propagating spin-0 
modes. Dynamical investigations using a simple expanding cosmological model found that the oscillation 
of the 0''' mode could account for an accelerating expansion similar to that presently observed. The model 
has been extended to include a 0^ mode and more recently cross parity couplings. We investigate the 
dynamics of this model in a situation which is simple, non-trivial, and yet may give physically interesting 
results that might be observable. We consider homogeneous cosmologies, more specifically, isotropic 
Bianchi class A models. We find an effective Lagrangian for our dynamical system, a system of first 
order equations, and present some typical dynamical evolution. 



1. Introduction 

All the known fundamental physical interactions can be fomiulated in a common framework: as local 
gauge theories. However the standard theory of gravity, Einstein's general relativity (GR), based on the 
spacetime metric, is a rather unnatural gauge theory. Physically (and geometrically) it is reasonable to 
consider gravity as a gauge theory of the local Poincare symmetry of Minkowski spacetime. A theory 
of gravity based on local spacetime geometry gauge symmetry, the quadratic Poincare gauge theory of 
gravity (PG, aka PGT) was worked out some time ago. 

We briefly sketch this theory, and how the search for good dynamical propagating modes led to 
focusing on the two scalar modes, and then how cosmological models were used to reveal the dynamics 
of the resulting model. 

There is no known fundamental reason why the gravitational coupling should respect parity. With 
this in mind, the general quadratic PG theory has recently had a major extension to include all possible 
odd parity couplings. The appropriate odd parity couplings have been incorporated, in particular, into the 
special case of the dynamically favored two scalar mode model. Here we show how a simple effective 
Lagrangian can reveal in a cosmological context the dynamics of this extended model. 



2. The Poincare gauge theory 

In the Poincare gauge theory of gravity |[1}|7], the two sets of local gauge potentials are, for "translations", 
the orthonormal co-frame -0°^ = e"jdx*, where the metric is 51 = — i?'' ® -d^ + dab'&"' ® i?^, and, for 
"rotations", the metric-compatible (Lorentz Lie-algebra valued) connection 1-forms F"'^ = T^"f^^idx\ 



The associated field strengths are the torsion and curvature 2-forms 

T" := di?" + A 1?'^ = ^r%^i9^ A 7?^, (1) 

^r"^ + r"^ A r^'^ = ^r'^^^.u^^' a r, (2) 

which satisfy the respective Bianchi identities: 

DT" = R",3 A DR'^fs = 0. (3) 

Turning from kinematics to dynamics, the PG Lagrangian density is generally taken to have the 
standard quadratic Yang-Mills form, which leads to quasi-linear second order equations for the gauge 
potentials. Qualitatively, 

^[i?, r] ~ K~^[A + curvature + torsion^] + fp^^curvature^, (4) 

where A is the cosmological constant, k = SirG/c^, and g^^ has the dimensions of action. The field 
equations obtained by variation w.r.t. the two gauge potentials have the respective general forms 

A + curvature + Dtorsion + torsion^ + curvature^ ~ source energy-momentum density, (5) 

torsion + Dcurvature source spin density. (6) 

From these two equations, with the aid of the Bianchi identities Q, one can obtain, respectively, the 
conservation of source energy-momentum and angular momentum statements. 

Earlier investigations considered only models with even parity terms; the models had in general 10 
dimensionless scalar coupling constants. The recent BHN investigation [Si systematically considered all 
the possible odd parity Lagrangian terms, introducing 7 new pseudoscalar coupling constants. Not all 
of these coupling constants are physically independent, since there are 3 topological invariants: the (odd 
parity) Nieh-Yan [9J identity 

A r„) = A + A 1?"^ (7) 

the (even parity) Euler 4-form i?"^ A R'^^rja/B^yS, and the (odd parity) Chern-Simons 4-form i?"^ A R^a- 
Early PG Investigations (especially [3, 10 1) of the linearized theory identified six possible dynamic 
connection modes; they carry spin-2^, spin-1^, spin-0^. A good dynamic mode should transport 
positive energy and should not propagate outside the forward null cone. The linearized investigations 
found that at most three modes can be simultaneously dynamic; all the acceptable cases were tabulated; 
many combinations of three modes are satisfactory to linear order. Complementing this, the Hamiltonian 
analysis revealed the related constraints [11]. Then detailed investigations of the Hamiltonian and 
propagation [12-15| concluded that effects due to nonlinearities in the constraints could be expected 
to render all of these cases physically unacceptable except for the two "scalar modes", carrying spin-0+ 
and spin-0^. 

In order to further investigate the dynamical possibilities of these PG scalar modes, Friedmann- 
Lemaitre-Robinson- Walker (FLRW) cosmological models were considered. Using a k = model it 
was found that the 0+ mode naturally couples to the acceleration of the universe and could account for 
the present day observations II16II171 ; this model was then extended to include the 0^ mode lITSl . 

After developing the general odd parity PG theory, in BHN [8| the two scalar torsion mode PG 
Lagrangian was extended to include the appropriate pseudoscalar coupling constants. 

The BHN Lagrangian jSl has the specific form 

1 ^ 1 
-2A + aoR + boX--Y^ a^^">r'^ + 3^2 T^P^ 



^BHN[^?,r] = ^ 



2 

n=l 



where R is the scalar curvature and X is the pseudoscalar curvature (specifically X/6 = i?[oi23] is the 
magnitude of the one component of the totally antisymmetric curvature), and T^j = r"a/t = '^'^^°'a^l, 
Pfi = ^Gfj.i'°'^T'^ai3 = are the torsion trace and axial vectors. 

In the aforementioned work the general field equations were worked out, and then specialized to 
find the most general 2-scalar mode PG FLRW cosmological model. Here we will consider manifestly 
isotropic Bianchi models. 

3. The PG scalar mode isotropic Bianchi I and IX cosmological models 

For earlier PG cosmological investigations see Minkevich and coworkers, e.g., |[T9l - l23l and Goenner & 
Miiller-Hoissen d; for recent work see ll8][T6HT8ll25ti28l. 

For the usual FLRW models, although they are actually homogeneous and isotropic, the 
representation is not manifestly so. Indeed they are merely manifestly isotropic-about-a-chosen-point. In 
contrast, the representation used here for the isotropic (class A) Bianchi I and IX models have the virtue 
of being manifestly homogeneous and manifestly isotropic. 

For homogeneous, isotropic Bianchi type I and IX (respectively equivalent to FLRW = and 
fc = +1) cosmological models the isotropic orthonormal coframe has the form 

:= dt, d"" := aa"", (9) 

where a = a{t) is the scale factor and depends on the (not needed here) spatial coordinates in such a 
way that 

da' = Ce'jk(7^ Aa^, (10) 

where = for Bianchi I and C = 1 for Bianchi IX, thus ("^ = k, the sign of the FLRW Riemannian 
spatial curvature. 

Because of isotropy, the only non-vanishing connection one-form coefficients are of the form 

^% = ^(t)a^ r\ = x{t)e\,a'^. (11) 

Here Cabc '■= ^abc] is the usual 3 dimensional Levi-Civita anti-symmetric symbol. 

From the definition of the curvature (|2j, one can now find all the nonvanishing curvature 2-form 
components: 

R\ = xdi A e V + [V' - x']ct'^ A (Tf, + xCe^ceV' A (12) 
= V'dt A ct'^ - xV'f^^ A e'^bcf^^ + VCeVc^^ A (T^ (13) 
Consequently, the scalar and pseudoscalar curvatures are, respectively, 

R = Q[a-^^ + a-\ij^ -[x-C? + C% (14) 

X = 6[a-'x + 2a-^^{x-C)]- (15) 
Because of isotropy, the only nonvanishing torsion tensor components are of the form 

T\o = u{t)dl T\, = -2x{t)e\c, (16) 

where u and x are the scalar and pseudoscalar torsion, respectively. From the definition of the torsion 
([T|l, one can find the relation between the torsion components and the gauge variables: 

n = a"^(d-V), X = a'^ix- 0- (17) 

Note that a~^ip = H — u, where H = a^^d is the Hubble function. 

From the isotropic assumption, the material energy momentum tensor must have the perfect fluid 
form. In our present work we focus on the late time behavior. Accordingly we assume that the fluid 
pressure can be neglected, so that the gravitating material behaves like dust with a density satisfying 
pa"^ = constant. Also, although we expect the spin density to play an important role in the early universe, 
it is reasonable to assume that the material spin density is negligible at late times. 



4. Effective Lagrangian 

The dynamical equations for these homogeneous cosmologies could be obtained by imposing the Bianchi 
symmetry on the general field equations found in [8| from their Lagrangian density. On the other 
hand dynamical equations can (following the same procedure used in 1 18 1) be obtained directly and 
independently from a classical mechanics type effective Lagrangian, which in this case can be simply 
obtained by restricting the BHN Lagrangian density to the Bianchi symmetry (this step is where the 
manifestly homogeneous representation plays an essential role). This procedure is known to successful 
give the correct dynamical equations for all Bianchi class A models (which includes our cases) in GR 
(see |[29l ). and it is conjectured to also work equally as well for the PG theory. Our calculations explicity 
verify this property for the isotropic Bianchi I and IX models. Indeed the equations we obtained in this 
way are equivalent to those found (at a later date) by BHN for their FLRW models (which are equivalent 
to our isotropic Bianchi models) by restricting to FLRW symmetry their general dynamical PG equations. 
This has proved to be a useful cross check. 
The ejfective Lagrangian is 



Leff = (2k) ^ {aoR + hoX - 2K)a^ + —i-a2u^ + 4.a-:>,x^ + 4.a2ux)a^ 

2k 



- —{wqR^ - W3X^ + fisRX) 



a . 



(18) 



It should be noted that the parameter restrictions 02 < 0, wq < 0, W3 > 0, and ^3 + Aw^wq < are 
necessary for the least action principle, which requires positive quadratic -kinetic-terms. 

In the following we often take for simplicity units such that k = 1 = g. These factors can be easily 
restored in the final results by noting that in the Lagrangian they occur in conjunction with certain PG 
parameters. Hence in the final results one need merely make the replacements 



{ao,a2,a3,bo,A,a2} ^ K ^{oq, 02, 03, 60, A, 0-2}, {w3,ws,^i3}^g ^{w3,Wfi, 1^3}. 



(19) 



The effective gravitational Lagrangian has the usual form of a sum of terms homogeneous in 
"velocities": L^s = Lq + Li + L2; the associated energy function is thus 



eff 



-W + -TT^X + - ^cff = L2- Lq 



die 



da 



'y -3(ao - -a2)u^ - 'iaoH'^ + 82:^(00 - 203) 

1 

+QuH{aQ - -02) + 6(5o + cF2)x{H - u) - 'ia^^ + A 
2 



^6 

24 

W3 



R^ - 12R \{H-uf -x'^ + ^ 



+_| [x' + 2AXx{H-u)] 



/f3 

24 



RX -6X UH -uy -x^ + -^} + 12Rx{H - u 



(20) 



Since L^q is time independent, the energy function is conserved. This constant energy has the form —a?p 
where p can be identified as the material energy density. This relation is consistent with the vanishing 
divergence of the source fluid energy-momentum tensor, which has been assumed to have the dust form 
with vanishing pressure. 



Making use of the formulas for the torsion and curvature components in terms of the gauge variables 
(I14I15I17I ). we now obtain the Euler-Lagrange equations. The dynamical equations are the if) equation: 

d dLes dLes 



2 4 



the X equation: 



--R+-X 

and the a equation: 



dt dip dip 
3(2ao - 02) - wqR - 



d dLes dLcs 



6(60 + ^2) - + 



X. 



dt dx 



dx 



6(60 + ^2) - Y^ + ^^sX 



u + 



6(ao - 203) - wqR - 



d dL 



cfT 



dL. 



cS 



dt da 



da ' 



(21) 



(22) 



— 3(a2U — 2cr2x) 



(ooi? + bfiX - 3A) + m{a2U - 2a2x) 

3 1 

+-(-a2M^ + 4032;^ + Aa2ux) - -^{wqR^ - w^X^ + fisRX) 



+ 



Oo 

2 
1 

24 



/•2 

-6(H-u? + 6x'^-6-^ 
a^ 



{2w^R + /i3X) 



-6(F - u) 



+ &)qx{H - u) 



(23) 



The above equations (I21I22I23I ) are 3 second order equations for the gauge potentials a, ip, x- 
However they can in an alternative way be used as part of a set of 6 first order equations along with 
the Hubble relation a = aH and the following two relations, obtained by taking the time derivatives of 
the torsion ([TT] ) and using the curvature definitions (I14ll5l i: 

X = -Hx 2x(H-u), (24) 

6 

H-ii = H(H -u)-{H -uf + x^ (25) 

6 

One advantage of such a reformulation is that the variables are now all observables. 

Our 6 first order dynamical equations and the energy constraint equation can now be put in the form 



a 
H 



x 



-"^R-^X 
2 4 



aH, 

1 
6a2 



d2R - 2a2X) - 2i?2 



Q2 - 4a3 2 _ C: , _P_ , 4A 



02 



3o2 3a2 



{aoR + d2X) - 3Hu + 



1 

3a2 

-^-{3H-2u)x, 

3(22 + WqR + ^X u + 

-6^2 + ^R - W3X 



4a3 2 , P , 4A 

^ + ^ + 

02 oa2 oa2 



^3 



u 



-6(T2 + —R - W3X 

1203 + weR + ^X 



(26) 
(27) 

(28) 

(29) 
(30) 
(31) 



3(-^a2 + 253)x2 + ^ 



A2' 



A 



Qa2 + W3X -^\x{H -u) + ^{weR^ - w^X^ + ^sRX) 



(H-u?-x'^ + -^ 



(32) 



where we have introduced certain modified parameters 02, 03 and (T2 with the definitions 



0,2 := 02 — 2ao, 



03 := 03 - y, a2 ■■= (72 + Oq. 



(33) 



The last two dynamical equations (I30I31I) could of course easily be resolved for R and X, but for some 
puposes the form given above is more convenient. 

It is remarkable that here in these effective Lagrange equations and the associated conserved energy 
we get (for this case with dust as the source) exactly the correct equations for our model — without 
including any explicit source coupling. 

From the above first order equations (following the method used in ifTSll ) we have found the associated 
linear equations and the late time normal modes. The details of these results will be reported in another 
work. In the following we will show some typical evolution. 



4.1. Numerical Demonstration 
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Figure 1. Evolution of the expansion factor a, the Hubble function, H, the second time derivative of the 
expansion factor, a, the scalar curvature and the pseudoscalar curvature, R and X, the torsion and the 
axial torsion function, u and x and the energy density, p. 



We present the results of a numerical evolution of our cosmological model. For all these calculations 
we take A = 0. We need to look into the scaling features of this model before we can obtain the sort 
of evolution results we seek on a cosmological scale. In terms of fundamental units we can scale the 
variables and the parameters as k = SvrG = 1. So the variables and the scaled parameters wq and 



^3 become dimensionless (from the Newtonian limit qq = 1). However, as we are interested in the 
cosmological scale to see changes of the age of our Universe, let us introduce a dimensionless constant 
To, which represents the magnitude of the Hubble time (Tq = Hq^ = 4.41504 x 10^^ seconds). With 
this scaling, all the field equations are kept unchanged while the period T — )■ TqT. 
We take the parameters as 

oo = 1, a2 = 0.83, ag = -0.35, wq = 1.1, W3 = 0.091, cj2 = 0.12 and /ig = -0.18. 

(34) 

The behavior of the 6 equations has been observed with several sets of initial values. We plot the typical 
case in (Fig. [U: first, evolution of the expansion factor a, second, the Hubble function, H, third, the 
second time derivative of the expansion factor, a, fourth, the scalar curvature and the pseudoscalar 
curvature, R and X, fifth, he torsion and the axial torsion function, u and x, sixth, the energy density, p. 

5. Concluding discussion 

We have been investigating the dynamics of the Poincare gauge theory of gravity. Recently, the model 
with two good propagating modes carrying spin 0^, 0^ (referred to as the scalar and pseudoscalar 
modes) has been extended to include pseudoscalar constants that couple the two different parity modes. 
It should be mentioned that here are many interesting physical and mathematical details related to torsion 
and parity and the type of models considered that we could not discuss here which are covered in 
references ||8][17l[T8l. Here we have considered the dynamics of this BHN model in the context of 
manifestly homogeneous and isotropic cosmological models. We found an effective Lagrangian and 
a system of first order equations. In these models, at late times the acceleration oscillates. It can be 
positive at the present time. It should be noted that the 0+ torsion does not directly couple to any known 
form of matter, but it does couple directly to the Hubble expansion, and thus can directly influence the 
acceleration of the universe. On the other hand, the 0~ couples directly to fundamental fermions; with the 
newly introduced pseudoscalar couphng constants it too can directly influence the cosmic acceleration. 
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